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Abstract 

Variational wave functions based on a Margenau-Brink cluster model with short 
range and state dependent correlations, and angular momentum projection are ob- 
tained for some nuclei with 12 < A < 16. The calculations have been carried 
out starting from the nucleon-nucleon interaction by using the Variational Monte 
Carlo method. The configuration used consists of three alpha clusters located at 
the apexes of an equilateral triangle, and an additional cluster, not necessarily of 
alpha type, forming a tetrahedron. This cluster is located at the top of its height. 
Short-range and state dependent correlations are included by means of a central 
Jastrow factor and a linear operatorial correlation factor respectively. Angular mo- 
mentum projection is performed by using the Peierls-Yoccoz operators. Optimal 
structures are obtained for all the nuclei studied. Some aspects of our methodology 
have been tested by comparing with previous calculations carried out without short 
range correlations. The binding energy, the root mean square radius, and the one- 
and two-body densities are reported. The effects of correlations on both the energy 
and the nucleon distribution are analyzed systematically. 

Key words: Nuclear structure; Cluster models; Variational Monte Carlo; N ^ Z; 
i>4 forces 

PACS: 21.60.-n; 21.60.Gx; 21.60.Ka; 27.20.+n 



* Corresponding author 

Email addresses: buendia@ugr.es (E. Buendia), galvez@ugr.es (F. J. Galvez), 
falsarua@uco.es (A. Sarsa). 



Preprint submitted to Elsevier Science 



8 February 2008 



1 Introduction 



The joint use of short-range dynamic correlations with model wave functions 
including relevant aspects of the nuclear structure constitutes the most com- 
monly used scheme to describe nuclear bound states with realistic or semi re- 
alistic interactions. Short range correlations are essential elements in the wave 
function because, as it is well known, any of the so-called realistic or semi- 
realistic parameterizations of the nuclear potential presents a strong short- 
range repulsive core. On the other hand, the formation of different kind of 
clusters in the nuclei can be understood as a collective movement of the nu- 
cleons governed by the medium and long range part of the nuclear potential. 
Therefore, for an accurate description of the nuclear states, it is convenient 
to consider both aspects in any variational approach to the nuclear bound 
states using this type of interactions. In principle, short range correlations are 
mainly governed by the nucleon-nucleon interaction while medium and long 
range effects depend on the particular nuclear state. However, and in a more 
careful approach, the final form of the short range correlations will depend 
on the model wave function giving rise to a non negligible dependence of the 
correlations on the nucleus. 

A direct way to include both short range and medium and long range cor- 
relations is by using Jastrow type correlation factors, but the calculation of 
the expectation values becomes very cumbersome, especially when state de- 
pendent correlations are included. There exist several methods to evaluate 
these expectation values as those based on cluster expansions [1,2], the Fermi- 
HiperNetted-Chain method [3,4] or statistical methods such as the Variational 
Monte Carlo [5,6]. The Coupled Cluster method allows to incorporate both 
type of correlations [7,8,9]. In this way it is possible to understand how the 
different correlation mechanisms are incorporated [10,11] 

Alpha cluster models, or cluster models in general, have been widely applied 
in microscopic descriptions of bound and scattering states of nuclear systems 
[12,13]. Variational wave functions built within this framework constitute an 
appropriate scheme for nuclei such as 8 Be and 12 C, that present a clear cluster 
structure. The use of wave functions including the possibility of formation of 
alpha cluster structures or any other kind of grouping of nucleons improves 
the description of these nuclei and their neighbours with respect to simple 
mean field approximations [14,15]. 

Multi cluster models have been used in microscopic calculations, i.e. with- 
out effective cluster-cluster interactions, based on the Generator Coordinate 
Method for some nuclei between A=12 and A=16 [16,17]. In these works a 
Volkov nucleon-nucleon potential was used [18]. Other results of microscopic 
multicluster calculations based on the stochastic variational method have been 



2 



reported [19,20] for some nuclei using the Minnesota potential. Neither of these 
potentials presents a strongly repulsive short range part and, therefore, short 
range correlations do not play a significant role. On the other hand, previous 
studies of alpha clustering based on nuclear potentials with a strongly repul- 
sive core have been mainly restricted to spin-isospin saturated nuclei [21,22]. 

The aim of this work is to study the ground state of some p-shell nuclei, 
A 7^ 4n, including clustering effects and short range and state dependent 
correlations, starting from t> 4 -type nucleon-nucleon interactions. The nucleon 
clustering is described in terms of model wave functions based on a general- 
ized Margenau-Brink model as in [16]. Short range correlations are included 
by means of a Jastrow factor and the dependence on the spin and isospin 
exchange channels is included by using a linear state dependent correlation 
factor. Angular momentum projection is carried out in order to obtain varia- 
tional wave functions that are eigenf unctions of the total angular momentum 
operator. The calculations are performed by means of the Variational Monte 
Carlo method. 

Here we extend a previous work [22] to the A ^ An case. This generalization 
is not straightforward because the angular momentum projection involves a 
spin mixing not present in spin and isospin saturated nuclei. In this paper 
we present an analytical reduction of the different expectation values for these 
nuclei, obtaining expressions suitable for the Variational Monte Carlo method. 
By using this scheme, the computing time is hardly increased with respect to 
the spin-isospin saturated case. We apply the method to the ground state of 
13 C, 14 C, 14 N, and 15 N. The results obtained are also valid for the mirror nuclei 
13 N, 14 and 15 because the electrostatic energy has been not considered in 
the minimization process. A systematic analysis of the effects of the different 
correlations mechanisms included in the wave functions on the total energy 
and on the contribution of the different channels is carried out. One and two 
body densities are reported and the effect of the correlations are discussed. 

The scheme of this work is as follows. In Section 2 the variational wave func- 
tion and the analytical reduction of the expectation values leading to a form 
appropriate for the Variational Monte Carlo method are detailed. In Section 3 
we report and discuss the main results here obtained. The conclusions of the 
present work can be found in Section 4. 



2 Wave function 

The variational trial wave function used in this work is 

*5 KM (1, 2, • • • , A) = Fj(l, A)F c (l, A)Q$ KM (1, ...,A) (1) 
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This structure has been used in previous studies of spin and isospin saturated 
nuclei [22,23]. It consists of a central Jastrow correlation factor Fj, a linear 
correlation factor Fc that can include state dependent correlations, and a 
model wave function ^ KM that is antisymmetric and has the proper values 
of the total angular momentum and parity. 

The Jastrow factor depends only on the distance between pair of nucleons 
Fj(l,...,A) = Y[f( rij ). (2) 

i<j 



The linear factor is defined as 

F £ (1,...,A) = X>(M) (3) 

i<j 



where the function g(i,j) depends on the radial and intrinsic degrees of free- 
dom of particles, This is the only part of the trial wave function where 
state dependent correlations are present explicitly. Here we employ the same 
parameterization for the correlation functions g(i,j) and f(r), used in previous 
works [21,22,23] that has shown to provide good results 

<K^) = i> (fc V«0P (fc Hu), (4) 

k=l 

where 

pW(i, j) = \ (l + n ■ tj) , = pW(i,j)pW(i,j). (5) 



This operatorial dependence of the correlation factor is the same as that of 
the nucleon-nucleon interactions considered in this work. The functions g^(r), 
k — 1, ..,4, and f(r) are parameterized as a linear combination of Gaussians 

M N 

9 {k \r) = E ^ e~ b - r \ f(r) = 1 + £ c n e~ d ^ . (6) 

m=l n=l 

The new aspects of treating A ^ 4n nuclei with respect to spin and isospin 
saturated ones are originated in the angular momentum projection. Therefore 
we shall focus here on the model part of the wave function and on the angular 
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momentum projection. The correlation factors are treated as in the spin and 
isospin saturated case. 

The model wave function used here is based on a generalization of the Margenau- 
Brink model. Instead of using only alpha-particle like nucleon clusters, more 
general groupings are allowed giving rise to a multicluster description [16,19]. 
Within the molecular viewpoint of the Margenau-Brink scheme, the model 
wave function is obtained starting from the following functions 

$c(l, 2, . . . ,A) = A {$i(rri, . . . ,x kl ) . . . ^ n {x kn _^ ...,x A )} (7) 



where C = {cfc})! =1 is a set of parameters that represent the centers of the 
clusters, and A is the corresponding antisymmetrizer. In this work the ar- 
rangement of the nucleons, shown in Fig. 1, consists of three a clusters and 
a fourth incomplete cluster that can be made of one, two or three nucleons 
depending on the nucleus under study. 




Fig. 1. Cluster description of the nuclei in terms of three alpha particles and a 
general s incomplete cluster with 1,2 or 3 nucleons 

For this configuration, the general form of the function given in Eq. (7) reduces 
to 



II $a m (£4m-3, • • • , %4r, 
m=l 



<f> s (x 13 ,...,x A )\ (8) 



where $ am stands for the wave function of an alpha particle centered at c, 
and $ s represents the incomplete cluster wave function centered at c s . 
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In this work the <J> Qm functions are taken to be Slater determinants built from 
harmonic oscillator single particle orbitals centered at c m 

= (y) 3 '* e-&^ (9) 

The oscillator parameter, /3, is the same for all of the alpha clusters. For the 
incomplete cluster wave function another Slater determinant centered at c s is 
employed also built from s-wave harmonic oscillator single particle orbitals. 
The oscillator parameter in this case is, in general, different to that for the a 
cluster wave function. The importance of using a different harmonic oscillator 
parameter will be discussed. With these choices for the cluster wave functions, 
the model wave function of the A nucleons is a Slater determinant. In gen- 
eral this function is not eigenfunction of parity or total angular momentum 
operators. 

The following linear combinations 

2, . . . , A) = $ e (l, 2, . . . , A) ± $_ e (l, 2, . . . , A) (10) 

have definite parity. Model wave functions with the total angular momentum 
of the state under study can be obtained from Eq. (10) by using the Peierls- 
Yoccoz projection operators [24] 

2 7 + 1 r 

*5 KM (1, ■ ■ ■ ,A) = j deV J M * K (Q)R(Q)^(l, . . . , A) (11) 

where R(O) is the rotation operator, T>j^ K {Q) is the rotation matrix and 9 
represents the Euler angles. The quantum number J gives the total angular 
momentum, K is its projection along the nuclear z axis and M the projection 
along the Z axis of the laboratory fixed frame. The projection within this 
scheme is carried out by rotating the intrinsic state and integrating over all 
angles weighted by the rotation matrix. 

The function • • . , A) in Eq. (8) is the generator function of the model 

wave functions. Note that we have removed the parametric dependence of the 
model wave function on the position of the centers, C, in order to simplify 
the notation. The distances between the clusters, R c and Rd, are determined 
variationally. 

The action of the rotation operator on the generator function is now described 
in detail. As we have mentioned before, this is the source of the new method- 
ological aspects originated by the fact that the nuclear states are not spin 
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and isospin saturated. We do not need to consider here the correlation fac- 
tors because they are rotationally invariant. The generator function is a Slater 
determinant. The action of the rotation operator on it leads to a linear com- 
bination of Slater determinants. If the Slater determinant is spin and isospin 
saturated this linear combination contains only one Slater determinant that 
also is spin and isospin saturated, containing the same single-particle orbitals. 
The only difference is that, after rotation, these orbitals depend on the rotated 
coordinates. This was exploited previously to study A = 4n nuclei [23]. When 
the nuclei are not spin or isospin saturated the rotation gives rise to a mixing 
of spin states. 

When the incomplete shell consists of one nucleon, as for example in the 
ground state of 13 C, the action of the rotation operator can be written as 
follows 



R(0)*a,^= E (12) 

where (5 stands for the spatial quantum numbers of the orbital of the incom- 
plete cluster, and sp and tp are the third component of spin and isospin, respec- 
tively. The over line indicates that the Slater determinant must be evaluated 
on the rotated coordinates. Therefore, and concerning to the spin dependence 
of the state, the effect of the rotation is to mix the two possible spin projec- 
tions of the orbital in the incomplete cluster. The weight of each component 
is given by the matrix element of the rotation matrix. 

When there are two extra nucleons the result of the rotation can be written 
as follows 



R ( Q )^c,p S0 t ,ps^ ~ 2 P 4 s i( e ) P ^( e )$cj3 Si t 5 ,/3 Sj t 7 

= 2 2 {\\sps 1 \S,sp + s, t ) (13) 
Si , Sj =±\ s=o,i zz 

(2 2~ SiS i\S, Si + Sj)^ , f /3+S7iS - +Si (Q)^ > c,/3s l t /3 ,/3s j t 7 

where {(3sptp) and (/?s 7 t 7 ) stand for the quantum numbers of the orbitals of the 
incomplete shell. Note that we have considered the same spatial dependence for 
both single particle orbitals. Therefore, if one is dealing with two extra protons 
( 14 0) or two extra neutrons ( 14 C) with the two possible spin orientations the 
term S = 1 vanishes. Only in the case of one proton and one neutron outside 
closed shell ( 14 N) both total spin components will contribute. 
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Finally, the case of three nucleons outside closed shell ( 15 N and 15 ) is a 
conjugate configuration to that of one nucleon outside closed shell and it is 
handled in the same way. 

The values allowed for J and K are governed by the symmetry group of the 
system, i.e. by the spatial positions of the centers of the clusters. For the 
nuclei here considered the group is C3„. The spin of the extra cluster must 
be also considered in determining the possible values of K. If Ms is the total 
spin third component the allowed K values are given by the selection rule 
\K — Ms\ = 3n, with n a positive integer [16], and, for any K, J > K and the 
parity is vr = (-1) J±S '. The energy grows with K , providing different rotational 
bands. In this work we are concerned only with the ground state, therefore 
we shall restrict ourselves to K — 1 for 14 N and K = for all the rest. For 
one and three extra nucleons Ms = 1/2 and the ground state is (l/2) + , and 
for two extra nucleons there are two possibilities; i) both nucleons are protons 
or neutrons M s = and the state is + , ii) one nucleon is a proton and the 
other a neutron Ms = 0, 1, and the 1 + ground state must be constructed with 
M s = 1 and K = 1. 

In order to compute the expectation value of the Hamiltonian in the projected 
wave function it is convenient to use the following expression [23,25] 

2 / + 1 r 

(^ M |H|^ M ) = J V J M * K (Q)(^\F c FjHFjF c K(Q)\<l>p(U) 



Let us focus on the spin-isospin configuration of the nuclear state. Note that, 
because of the rotational invariance property of the Hamiltonian, only the ket 
is rotated remaining the bra on its original configuration. This is important 
because it determines the configurations that gives non zero contribution to 
the integral when projected onto the bra. The action of the rotation operator 
is to produce a linear combination of configurations containing the original 
one. One needs to analyze all of them to determine if, after the action of the 
spin-isospin operators of Fc and the Hamiltonian, the original configuration is 
obtained. As a result, only the original configuration appearing after rotation 
contributes with both central and state dependent correlation factors, except 
except for incomplete clusters made of one proton and one neutron with S = 0, 
that we have not studied here, for which two of the configuration appearing 
after rotation give non zero contribution. Note that the weight factor must be 
included when doing the integral in all of the cases. The treatment of state 
dependent correlations in terms of the intermediate states is not modified with 
respect to the case of spin and isospin saturated nuclei [26,27]. 
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Table 1 

Binding energy and root mean square radius, (r 2 ) 1 / 2 , for different nuclear states 
calculated in this work (mc) as compared with the results of Dufour and Descou- 
vemont (dd) [16]. Both calculations have been performed by using the Volkov V7 
interaction [18] and the same variational wave function without correlations. The 
inverse of the oscillator parameter, and the distances between the clusters, 

R c and Rj, are also included. The energies are in MeV, and (r 2 ) 1 / 2 , /3 _1 , and R c 
and Rd, in fm. The statistical error in the Monte Carlo calculation is indicated in 
parentheses. The Coulomb energy has been included in the total energy. 



A X(K, J w ) 


r 1 


Rc, Rd 


Emc 


Edd 


/_2\l/2 
V /mc 


c^ 2 


12 C(0,0+) 


1.38 


2.65 


86.49(4) 


86.7 


2.31(7) 


2.31 


12 C(3,3-) 


1.38 


3.14 


76.41(4) 


76.5 


2.49(9) 


2.49 


13 c(U") 


1.39 


2.29,2.114 


88.99(7) 


89.6 


2.25(9) 


2.25 


14 C(0,0+) 


1.39 


2.26,2.057 


102.26(6) 


102.5 


2.26(7) 


2.26 


15 N(I,I-) 


1.35 


1.84,1.887 


119.37(7) 


121.9 


2.15(11) 


2.15 


16 O(0,0+) 


1.34 


1.49,2.409 


147.83(5) 


148.0 


2.18(3) 


2.18 


16 0(3,3-) 


1.37 


2.24,1.958 


129.46(10) 


129.8 


2.27(10) 


2.26 



3 Results 



First we will test the new methodological aspects implemented in this work 
by comparing with the results of Dufour and Descouvemont [16] obtained by 
using a different computational scheme. We will employ for the test both the 
same nucleon-nucleon interaction (the Volkov V7 potential), and the same 
wave function as in [16]. It is worth to point out that the correlation factor is 
not needed because the interaction does not present a strongly repulsive core. 

In Table 1 we show for the ground state and some excited states of the nuclei 
studied in this work the binding energy and the root mean square radius, 
(r 2 ) 1//2 . As can be seen from the table, both set of results are in a very good 
agreement. The spin-orbit interaction is not included in our work and therefore 
one can not compare directly the results for nuclei with an odd number of 
nucleons. For these nuclei we have compared with the average value of the 
states l/2~ and 3/2~ of [16]. This average gives a value that it is very close to 
the Monte Carlo result of this work, specially for 13 C where the spin-orbit 
splitting is smaller than in 15 N. From this test it can be concluded that, 
for A ^ An, the angular momentum projection scheme of this work provide 
reliable results. 

The ground state of these nuclei has been studied in this work by using a semi 
realistic potential. We have used the modified Afnan-Tang nuclear potential 
MS3 [28,29]. This is a v 4 type interaction with a strongly repulsive core. It 



9 



gives meaningless results when used with non correlated trial wave functions. 
Thus, in order to analyze the effects of nuclear correlations with respect to 
the non-correlated case, it is more convenient to use an interaction with a less 
repulsive short range part as the Brink-Boeker BB1 force [30]. 

The ground state energy and the root mean square radius (r 2 ) 1 / 2 for different 
nuclei calculated from a number of trial wave functions by using the BB1 
and the MS3 interactions are reported in Tables 2 and 3, respectively. The 
optimal parameters of the trial wave functions are also shown. The notation 
is as follows: MB stands for a non correlated trial wave function, JL includes 
central Jastrow and linear state independent correlations and JLO is a wave 
function with central Jastrow and state dependent linear correlations. In the 
JLO approach we have used the same non-linear parameters as in JL, i.e. the 
variational freedom is restricted only to the linear parameters of the different 
operatorial channels. This scheme has shown to work properly for spin and 
isospin saturated nuclei [22,23] in such a way that the loss of energy due to 
this partial optimization was very small. This is convenient because when state 
dependent correlations are included, two things happens; first the calculation 
becomes slower, and second, the statistical error increases. Therefore it is very 
convenient, from a computational point of view, that the non-linear parameters 
can be well determined by means of a state independent optimization. Note 
that the linear parameters are computed by solving a generalized eigenvalue 
problem and then only a long run is required to fix them. The expectation value 
of the Coulomb energy E c , not included in the total binding energy, is reported 
separately. For the results shown in this work we have used 2 s x 10 5 (2 6 x 10 5 ) 
moves per-nucleon with state independent (state dependent) correlated wave 
functions. 

The wave functions used in this work includes two different oscillator param- 
eters, one for the complete clusters and another for the incomplete one. This 
gives rise to an improvement in the energy of about 3 or 4 MeV when the 
incomplete cluster is made of one or two nucleons. The improvement is notice- 
ably reduced if the incomplete cluster contains three nucleons. The smaller 
value for the oscillator parameter of the incomplete cluster is due to the fact 
that the nucleons are more localized in the alpha particle cluster than in the 
incomplete cluster. In general we have obtained oscillator parameters that 
vary between those of 12 C and 16 0. 

With respect to the optimum parameters of the inter-cluster distances, we 
have obtained that the distance between the centers of the complete clusters 
is bigger than the distance between the incomplete cluster and an alpha- 
particle cluster. The total energy is not very sensitive to variations of the 
inter-cluster distances in the neighbourhood of the equilibrium values. We 
have indicated such situation by giving these distances with only one decimal 
digit. Finally and, as it could be expected, when moving from A = 12 to 
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Table 2 

Ground state energies calculated by using different trial wave functions without 
correlations (MB), with state independent correlations (JL) and with linear state 
dependent correlations (JLO) for the BB1 Brink-Boeker potential. Energies are in 
MeV, (r 2 ) 1 / 2 in fm, /3i,/?2 in fm -1 and R c ,Rd, in fm. The statistical error is shown 
in parentheses. The Coulomb energy is not included in the total energy. 



A(ii,J ) 


WTP 


ft-. Ft 
Pi 7 P 


2 


E> 




TP 
Hi 




/ r 2\l/2 
V ) 




1V1H) 


n 7n 

u. / u 




o.^± 




8D D1 ( A\ 


7 1 Q7( 1 "\ 
t .iy 1 11 j 


9 K'itA} 

Z.UO^ J 




TT 
JL 


n 70 
U. 1 1 




O.O 




— 112.00(4 ) 


( .41 f (1J 


z.0o( 1 j 




TT O 

J 1 j v 7 


n 79 

u. / Z 




O.O 






/ .oy / yij 


9 ^1(7) 

Z.OO^ / J 




T\/TO 
IVlO 


U.Oo, 


u.oy 


O.O, 


O.U 


78 IQffW 
— 1 O.Zy (O ) 


1 .UO 1 \ J-J 


9 71 (Q\ 




TT 

JL 


u. r z, 


U.54 


O A 

3.4, 


3.0 


— llZ.OO( 1 ) 


^.ooo(l J 


z.o3(o J 




TT n 


n 79 
u. / z, 




0.41:, 


O.U 


— 1 iy. o(z j 


1 .U10(Z J 


Z.OZ^lO J 




A/TO 


u.oy, 


U.OD 


9 

o.z, 


9 n 
z.o 


— oO.OO(0 J 


( .oOo(l ) 


9 fi/1 




TT 

JL 


U. 74, 


U.OO 


3.1, 


c 
z.o 


— lzz.93(8J 


/.oo6(l J 


z.4^(oj 




TT n 


n 7zL 
U. / 4L, 


U.Oo 


0.1, 


9 8 
z.o 


— 101. / O^IO J 




9 ARl'S^ 
Z.^U(o J 




MB 


0.68, 


0.57 


3.2, 


2.8 


-85.09(6) 


9.849(1) 


2.65(8) 


IN (1, 1 J 


TT 
JL 


U. / 1, 


U.O t 


o.U, 


k 
z.o 


— Izl.oo^ i ) 




2.47(7) 




TT O 


71 


u.o / 


O.U, 


Z.O 


— 101 .O ^Z J 


1 n SRi ^9^ 

lU.OOl ^Z ) 


Z.^rO^lU J 




A/TO 


U.OO, 


U.OD 


o.U, 


9 ^ 

Z.O 


— u i .uy ^iu j 




z.oo(y j 


1 ^ TV T / 1 1 — \ 


JL 


0.74, 


0.63 


2.7, 


2.4 


-139.55(10) 


10.821(1) 


2.39(9) 




JLO 


0.74, 


0.63 


2.7, 


2.4 


-152.0(4) 


10.837(5) 


2.38(18) 




MB 


0.66 




2.9, 


2.4 


-118.70(5) 


13.470(1) 


2.60(3) 


16 O(0,0+) C3 „ 


JL 
JLO 


0.76 
0.76 




2.8, 
2.8, 


2.4 
2.4 


-166.92(6) 
-179.46(10) 


14.516(1) 
14.515(2) 


2.36(3) 
2.35(5) 




MB 


0.67 




2.8 




-118.52(5) 


13.456(1) 


2.60(3) 


16 O(0,0+) t 


JL 


0.74 




2.6 




-166.66(6) 


14.446(2) 


2.37(4) 




JLO 


0.74 




2.6 




-180.61(8) 


14.552(2) 


2.35(5) 



A = 15 the optimal values of the variational parameters tend to those of O. 
This is the case for all of the interactions and wave functions analyzed in this 
work. It is remarkable that the ground state energy of 16 obtained with the 
C 3v symmetry is practically the same as the one obtained with a tetrahedral 
symmetry. 

In general, the effect of the correlations is to reduce the average size of the 
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Table 3 

Ground state energies calculated by using different trial wave functions without 
correlations (MB), with state independent correlations (JL) and with linear state 
dependent correlations (JLO) for the modified Afnan-Tang MS3 potential. Energies 
are in MeV, (r 2 ) 1//2 in fm Pi, 02 in fm -1 and R c ,Rd, in fm. The statistical error is 
shown in parentheses. The Coulomb energy is not included in the total energy. 



A X(K, J w ) 


WF 


01,02 


Rc, Rd 


E 


E c 


( r 2)l/2 


12 C(0,0+) 


JL 


0.70 


3.5 


-74.54(5) 


7.571(1) 


2.48(4) 




JLO 


0.70 


3.5 


-87.2(4) 


7.440(2) 


2.49(15) 


13 c(iO 

v 2 ' 2 / 


JL 


0.70, 0.46 


3.3, 3.1 


-73.37(10) 


7.833(1) 


2.47(8) 




JLO 


0.70, 0.46 


3.3, 3.1 


-88.6(6) 


7.864(1) 


2.44(13) 


14 C(0,0+) 


JL 


0.69, 0.48 


3.4, 3.0 


-77.52(7) 


7.840(1) 


2.50(5) 




JLO 


0.69, 0.48 


3.4, 3.0 


-94.6(3) 


7.840(1) 


2.44(10) 


14 N(1,1+) 


JL 


0.69, 0.54 


3.2, 2.8 


-81.95(9) 


10.699(1) 


2.42(7) 




JLO 


0.69, 0.54 


3.3, 3.8 


-99.3(4) 


10.865(3) 


2.37(10) 


15 N(I,I-) 


JL 


0.67, 0.54 


3.2,2.8 


-91.77(12) 


10.701(1) 


2.45(9) 




JLO 


0.67, 0.54 


3.2, 2.8 


-112.6(6) 


10.878(3) 


2.39(15) 


16 O(0,0+) 


JL 


0.71 


2.7 


-114.46(7) 


14.827(1) 


2.32(3) 




JLO 


0.71 


2.7 


-135.6(3) 


15.036(2) 


2.27(7) 



nucleus. Therefore, the optimum values in the model wave function will de- 
pend on the presence, or not, of the correlation factor. The modification with 
respect to the non-correlated wave function is roughly proportional in all of 
the parameters in such a way that nucleon correlations give rise to an isotropic 
contraction of the nucleus. 

It is interesting to point out the importance of correlations in the binding 
energy of 12 C and 14 C as compared with 13 C and 14 N, respectively. With both 
interactions, 12 C is more bounded than 13 C with central correlations, but state 
dependent correlations reverse this situation, obtaining a difference of 1 and 
2 MeV with the MS3 and BB1 interaction, respectively. The behaviour of the 
nuclear binding energy of 14 C and 14 N is different with both potentials. With 
the BB1 interaction, and without correlations, 14 C is slightly more bounded 
than 14 N. The difference in their binding energy decreases with the use of cen- 
tral correlations and is zero with state dependent correlations. However, with 
the MS3 potential, 14 N is 4.5 MeV more bounded than 14 C with central and 
state dependent correlations. The reason of this different behaviour lies in the 
contribution of the Bartlett and Heisenberg channels of the MS3 interaction, 
that are null in the BB1 potential. Finally it is also worth mentioning here 
that we have obtained a negligible effect of the state dependent correlations on 
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Table 4 

Increase in the binding energy per number of nucleon pairs due to the inclusion 
of different correlation factors for the nuclei studied in this work. In parentheses is 
indicated the nuclear interaction. The increment is in MeV per number of nucleon 
pairs. The error is in the last figure. 



A X(K,J*) A JL _ MB (BB1) A JLO -jl(BB1) A JLO -jl(MS3) 



12 C(0,0+) 


-0.49 


-0.08 


-0.19 


i3 c(in 


-0.44 


-0.09 


-0.19 


14 C(0,0+) 


-0.40 


-0.10 


-0.19 


14 N(1,1+) 


-0.40 


-0.11 


-0.19 


15 N(I,I-) 


-0.40 


-0.12 


-0.19 


16 O(0,0+)(C 3 „) 


-0.40 


-0.12 


-0.18 



the Coulomb energy, which depends basically on the parameters of the model 
wave function. 

The correlations increase the binding energy by a quantity which grows with 
the number of nucleons A. In order to get a deeper insight into the coupling 
between correlations and the particular nucleus we report in Table 4 the in- 
crement in energy per number of pairs of nucleons. For example the increase 
in the binding energy per nucleon pair when state independent correlations 
are included with respect to the uncorrelated model is given by 

2 

Ajl-MB = ~T7~ A 7T (-^JL — -^Mb) 



where (£mb) is the energy in the JL (MB) model. The quantity Ajlo-jl 
is denned in a similar way. As it can be seen, the increment per number of 
pairs is roughly constant for all of the nuclei considered, specially Ajlo-jlj 
that accounts for the effect of state dependent correlations. The increment due 
to state dependent correlations in the MS3 potential is practically twice the 
increment in the BB1 case. 

A more detailed analysis of the effect of the state dependent correlations on 
the energy can be done by looking at the contribution of the kinetic energy 
and of the different channels of the potential energy. In Fig. 2 we plot the 
differences between these quantities calculated with the JL and JLO wave 
functions for both the BB1 and the MS3 interactions. Both the kinetic energy 
and the energy of the Wigner channel rise with state dependent correlations 
for both potentials. This increase is more important for the kinetic energy 
with the MS3 potential than with the BB1 one, whereas the opposite holds 
for the energy of the Wigner energy. For the BB1 potential, the Majorana 
channel is the responsible for the decrease in the ground state energy when 
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state dependent correlations are considered. For the MS3 interaction, the effect 
on the Majorana channel is practically canceled with that on the kinetic and 
Wigner energies, and the Bartlett and Heisenberg channels make the nuclei 
more bounded. The contribution of these two channels is very close and is 
nearly independent of the nucleus considered. 





11 a x a a 11 a x a a 

.■i .-i 

Fig. 2. Increase in the total energy, the expectation values of the kinetic energy 
and the different channels of the interacting potential when state dependent are 
included with respect to the JL approximation. In the left hand panel we plot the 
results for the BB1 potential and in the right hand one for the MS3 potential. The 
lines are for guiding the eyes 

The one- and two- body densities give the spatial distribution of the nucleons 
in the nuclei. Here we have calculated these densities to analyze the effect of the 
different correlation mechanisms introduced in the variational wave functions. 
In Fig. 3 we show the one body nuclear density calculated with the JL wave 
function for all of the nuclei here studied and the two interactions considered. 
As it can be seen, the qualitative behaviour is similar for both potentials, with 
a higher value of the maximum as the number of nucleons increases. It is also 
worth pointing out that as A increases, the density tends to that of 16 0. It is 
for this nuclei and the MS3 interaction where this density is more separated 
from the others. 

The effect of the state dependent correlations on the one-body density for 
these nuclei is studied in Fig. 4, where we plot the difference between the 
single particle density obtained with the JL and the JLO wave functions. The 
first noticeable fact is that the general behaviour is different for the two inter- 
actions used here. Thus at short distances state dependent correlations tend to 
increase the density with the BB1 interaction and the opposite happens with 
the MS3 potential, except for 12 C for which a negative region at short distances 
appears. In addition, for the BB1 potential, the effect of the operatorial cor- 
relations is roughly independent of the nucleus while for the MS3 potential 
effects of the operatorial correlations show a more accused dependence on the 
nucleus. 

The effects of correlations are more important on the two body density than 
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Fig. 3. One body density for all the the nuclei studied in this work calculated 
with the JL wave function. In the left hand panel we plot the results for the BB1 
potential and in the right hand one for the MS3 potential. 




r lltajl r lltajl 



Fig. 4. Effect of the state dependent correlations on the one body density for the 
different nuclei considered in this work. In the left hand panel we plot the results 
for the BB1 potential and in the right hand one for the MS3 potential. 

in the one body density. In Fig. 5 we plot the two body density obtained from 
the state independent correlated wave function JL for all of the nuclei studied 
and the two interactions considered in this work. The behaviour of this density 
is very similar for both potentials, although the effect of the nuclear core is 
much more important in the MS3 potential. The main difference is that with 
the MS3 interaction shorter distances are favoured with respect to the BB1 
potential. At distances between 2 and 3.5 fm the differences among the nuclei 
considered are more important, with bigger values as the number of nucleons 
increases from 12 C to 16 0. This can be understood as a progressive filling of 
the incomplete cluster that gives rise to a larger number of particles at these 
intermediate distances. 

Finally, the effect of including state dependent correlations on this density is 
studied in Fig. 6 where we plot the difference between the two-body density 
calculated from the JL and JLO wave functions. As it was the case for the 
one body density, the effect of state dependent correlations is roughly inde- 



15 



d 1 1 ■ £ ■ ■ d 1 1 ■ £ 

r iWl r iWl 



Fig. 5. Two body density for all the the nuclei studied in this work calculated 
with the JL wave function. In the left hand panel we plot the results for the BB1 
potential and in the right hand one for the MS3 potential. 

pendent of the nucleus when the BB1 potential is used and a more accused 
dependence is observed for the MS3 interaction. For both potentials, state 
dependent correlations bring together nucleons with respect to the JL case. 




Fig. 6. Effect of the state dependent correlations on the two body density for the 
different nuclei considered in this work. In the left hand panel we plot the results 
for the BB1 potential and in the right hand one for the MS3 potential. 

4 Conclusions 

Variational Monte Carlo calculations for p-shell, A ^ An, nuclei starting from 
the nucleon-nucleon interaction have been presented. The ground state energy 
and the one and two body densities have been calculated. The variational wave 
function consists of three factors: a central Jastrow term, a spin-isospin de- 
pendent linear term and a model wave function. The model wave function is 
based on a cluster model allowing for the formation of different kind of nucleon 
clusters with centers at fixed positions. The Peierls-Yoccoz projection opera- 
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tors have been used in order to obtain trial wave functions with the proper 
values of the angular momentum. This work extend previous ones carried out 
for spin and isospin saturated nuclei. 

The present scheme has shown to be appropriate for describing two important 
and complementary aspects of the nuclear dynamics as the short range cor- 
relations and the formation of nucleon clusters. The former is induced by the 
short range repulsive part of the nuclear potential while the later is a collective 
effect due to the medium and long range part of the interaction. 

In this work, an analytical reduction of the expectation values for A ^ An 
nuclei is presented. The use of the Peierls-Yoccoz projection operators intro- 
duces new features when the nuclei are not spin and isospin saturated. Here we 
obtain a final form of the expectation values which is specially suited for the 
Variational Monte Carlo calculation. This is done for both state independent 
and state dependent correlation factors. As a result the different expectation 
values can be computed with no significant extra computational cost with 
respect to the case of spin and isospin saturated nuclei. 

The scheme is applied to several nuclei with 12 < A < 16. The methodol- 
ogy has been first tested against previous works using a completely different 
scheme of calculation. Then results obtained by using two different nucleon- 
nucleon potentials including a repulsive core at short distances and state- 
dependent interaction channels have been reported. The binding energies and 
the root mean square radius along with the optimal parameters of the wave 
functions are shown for the different nuclei and states considered here. The ef- 
fect of the different correlation mechanisms included in the trial wave function 
on the energy and on the equilibrium geometries is discussed. The importance 
of using different oscillator parameters for the different kind of nucleon clusters 
is shown. The effect of the correlations on the different interaction channels 
is analyzed in terms of the number of nucleons. Finally one and two body 
densities obtained for the nuclei here studied with several approximations of 
the wave functions are reported and discussed. 
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